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Abstract 



The data arrangement problem on regular trees (DAPT) consists in assigning 
the vertices of a given graph G to the leaves of a d-regular tree T such that the sum 
(~| ' of the pairwise distances of all pairs of leaves in T which correspond to edges of G 

is minimised. Luczak and Noble [B] have shown that this problem is A^P-hard for 
every fixed d >2. 

In this paper we propose construction and local search heuristics for the DAPT 
and introduce a lower bound for this problem. The analysis of the performance of 
the heuristics is based on two considerations: a) the quality of the solutions produced 
^ ' by the heuristics as compared to the respective lower bounds b) for a special class 

of instances with known optimal solution we evaluate the gap between the optimal 
Qs ' value of the objective function and the objective function value attained by the 

l/^ , heuristic solution, respectively. 

\l ■ Keywords. Combinatorial optimisation; data arrangement problem; regular trees; 

^^ ' heuristics. 
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^ ; 1 Introduction 



Given an undirected graph G = (y{G),E{G)) with |1^(G)| = n, an undirected graph H = 
{V{H),E{H)) with \V{H)\ > n and some subset B of the vertex set of H, B C V{H), 
with \B\ > n, the generic graph embedding problem{GF,P) consists of finding an injective 
embedding of the vertices of G into the vertices in B such that some prespecified objective 
function is minimised. Throughout this paper we will call G the guest graph and H the 
host graph. A commonly used objective function maps an embedding (p: V{G) -^ B to 

E d{<P{i)cP{j)) , 
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where d{x, y) denotes the length of the shortest path between x and y in H. The host 
graph H may be a weighted or a non-weighted graph; in the second cases the path lengths 
coincide with the respective number of edges. Given a non-negative number A E M, the 
decision version of the GEP asks whether there is an injective embedding 0: V{G) —?■ B 
such that the objective function does not exceed A. 

Different versions of GEP have been studied in the literature; the linear arrangement 
problem, where the guest graph is a one dimensional equidistant grid with n vertices, 
see [21 [5l |9], is probably the most prominent among them. A number of other classical 
and well known combinatorial optimisation problems can be seen as special cases of 
the GEP, as e.g. the Hamiltonian cycle problem, the Hamiltonian path problem and the 
graph isomorphism problem (see e.g. [l] for a more detailed discussion of the relationship 
between these problems). 

This paper deals with the version of the GEP where the guest graph G has n vertices, 
the host graph H is a, complete d-regular tree of height [log^ n] and the set B consists 
of the leaves of H. From now on we we will denote the host graph by T. The height of 
T as specified above guarantees that the number |i?| of leaves fulfills \B\> n and that 
the number of the predecessors of the leaves in T is smaller than n. Thus [log^^ n\ is the 
smallest height of a d-regular tree which is able to accommodate an injective embedding 
of the vertices of the guest graph on its leaves. This problem is originally motivated by 
real problems in communication systems and was first posed by Luczak and Noble [6]. 
We will call this version of the GEP the data arrangement problem on regular trees 
(DAPT). Luczak and Noble [6| have shown that the DAPT is A^P-hard for every fixed 
d > 2. The question about the computational complexity of the DAPT in the case 
where the guest graph is a tree, posed by Luczak and Noble in [6], is still open. In this 
perspective the development of heuristic approaches to efficiently find good solutions to 
DAPT is a natural task. There are plenty of heuristics for different versions of the GEP 
in the literature, especially for the linear arrangement problem, see e.g. the papers by 
Petit [TJ E] for nice and comprehensive reviews. However, to our knowledge there are 
no specific heuristic approaches to solve the DAPT and no benchmark instances have 
been developed for this problem yet. In this paper we make a first step in this direction 
and propose construction and local search approaches as well as a lower bound for the 
DAPT, much in the spirit of [71 E] which deal with the linear arrangement problem. In 
order to evaluate the performance of the proposed heuristics we generate a number of 
families of test instances some of them being polynomially solvable or having a known 
optimal objective function value. 

The paper is organised as follows. Section [2] discusses some general properties of 
the problem and introduces the notation used throughout the paper. In Section [3] we 
derive a lower bound for optimal objective function value to be used in the evaluation 
of the performance of solution heuristics. Section [J] introduces the proposed heuristics. 
Sections [5l [6] and [7] discuss the test instances, the numerical results and some conclusions 
and outlook, respectively. 



2 Notations and general properties of the DAPT 

Consider a guest graph G = {V, E) with n vertices, \V\ = n, and a host graph T which 
is a d-regular tree of height h, h := flog^ n] . Let B be the set of leaves of T. Notice that 
due to the above choice of h we get the following upper bound for the number b = \B\ 
of leaves: 



b:= \B\ =d^ = d!''^d<nd. 



(1) 



Definition 1 An arrangement is an injective mapping 4> : V —?■ B. The data arrange- 
ment problem on regular trees (DAPT) asks for an arrangement that minimises the 
objective value OV{G, d, (p) 



OV{G,d,c^):= Yl dTiHu),Hv)) 



(2) 



(u,v)£E 



where dT{4>iu),(p{v)) denotes the length of the (f){u)-(j){v)-path in the d-regular tree T. 
Such an arrangement is called an optimal arrangement. An instance of the DAPT is 
fully determined by the guest graph and the parameter d of the regular tree T which serves 
as host graph. Such an instance of the problem will be denoted by DAPT{G, d) . 

Figure m shows a guest graph G with vertices {vi.,V2.,v^.,Vi.,v^} and Figure [2] shows a 
3-regular tree of height 2 = [log3 5] as a host graph together with a minimum arrange- 
ment. The numbers in the leaves of T denote the vertex indices mapped to the leaves, 
respectively, 




Figure 1: A guest graph. 




Figure 2: An optimal arrangement which objective value is 20. 



Notations. From now on let the set of vertices of the guest graph G be given as V{G) = 
{vi, . . . Vn} and let m := \E\ be its number of edges. We denote the set of neighbours of 
any vertex v by T{v). We denote by h{T) the height of a ((i-)regular tree T. A basic 
subtree T' of the d-regular tree T is a d-regular subtree of T with h{T') = h{T) — 1 rooted 
at some son of the root of T. For every d,h gN, 2 < d < n, the leaves of a d-regular tree 
of height h are denoted by bi,b2, ■ ■ ■ ,b^h such that bu_i\^h 



-1+1) ^(j- 



l)d'>-i+25 • • • '"«•( 



are the leaves of the i-th basic subtree. This order of the leaves is called the canonical 
order. If the leaves are labelled according to the canonical ordering then the pairwise 
distances between the leaves of a d-regular tree are given by a simple formula. 

Observation 2 Let T be a d-regular tree of height h and let its leaves be labelled ac- 
cording to the canonical ordering. The distances between the leaves in T are given as 
dT{bt,bj) = 21, where 



I :-- 



mm 



ke {1,2,. ..,h}: 



t-1 

[ d^ \ 


= 


i-1 



for all leaves bf, bj of T with t,j G {1, 2, . . . ,d^}. 



First let us observe that for aU t,j € {1, 2, . . . , d^}, [^ 



Proof. 

[^\ = [jwr\JoTallle{l,2,...,h-l}. 
We prove the claim by induction on h 



[^-jr] implies 



If /i = 1 then T has d leaves labelled by 



1,2, ... ,d, their pairwise distances are all equal to 2 and L^^J = L'^^J = 0) so the claim 
holds. Assume that the claim holds for regular trees of height up to h — 1. Consider now 
a tree of height h with leaves labelled by bi,b2, . . . , b^h in the canonical ordering and let 
bt, bj be two leaves of it. 



Let t = {it — l)d^ ^ + Tt and j = {ij — l)d^ ^ + rj with it,ij G {1, 2, . . . , d} and 
rt,rj G {l,2,...,d^-i}. Clearly \j^] = it - 1 and [^^1 = ij - 1. Thus, if I := 

min{/c G {1,2, ... ,h}: L^^J = L'^d^Ji — ^' ^^ ^^^ '^^^V ^^ ^t t^ ^i' ^^ equivalently, bt, bj 
are leaves of different basic subtrees of T. For leaves of different basic subtrees we have 
dT(bi,bj) = 2h and hence, the claim holds in this case. 



Otherwise / := min{A; G {1,2, . . . ,h}: [^J = [^J } <h-l which implies [/ ^ 



[4jt7tJ • Thus bt and bj are leaves of the same basic subtree of T. Let this be the r-th basic 
subtree T^ of T with leaves ^(r-i)d''-i+s with s = 1,2,..., d^~^. In the canonical ordering 
in Tr these leaves would be labelled by bg, for s = 1,2,..., d^^^. Let t = [r — l)d^^^ + st 
and J = (r — l)d 

hence dr^ (bst , bs. 



+Sj for St, Sj (^ {1,2, ... ,d }. T^ is d-regular tree of height h—1 and 
= 21 holds, where I := mm{k e {1,2, . . . ,h - 1}: [^\ = [^J }, 

dxipt, bj) and 



according to our inductive assumption. Finally notice that dTribstibs/, 



mm <; A; G {1,2, 
mm <, k £ {1,2,. 



1}: 
1}: 



st-1 

d^ 
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Sj — 1 
d^ 


} 




(r - l)d^-^ + St - 1 


= 




d^ 







l)d' 



h-l 



+ Sj 



d^" 



hold. This completes the proof. D 

Definition 3 For a given arrangement (p let Bu = {0(1), • • • , 4'i''^)} be called the set of 
used leaves. // Bu = {ftj, . . . &j+„_i} holds for some 1 < i < b — n + 1, cf) is called a 
contiguous arrangement. 

Let us notice that not every instance of the DAPT possesses necessarily a contiguous 
optimal arrangement as illustrated by the following example. 

Example 4 A DAPT instance which does not possess any contiguous optimal arrange- 
ment. 

The guest graph G with 12 nodes is represented in Figure Consider d = A. The 
optimal arrangement (j) represented in Figure [^ is not contiguous. In both pictures we 
identify the vertices with their indices, thus we write i instead of Vi, i = 1,2, ... , 12 for 
simplicity. The optimal value OV{G, 4, (p) equals 28 and can be written as OV{G, 4, (p) = 
4:*a{(p)+2{m—a{(p)), where m = 11 is the number of edges of the guest graph and a{<j)) = 3 
is the number of edges of G with end-vertices mapped by cp into different basic subtrees 
ofT. 

We show now that for every contiguous arrangement tp, a{'ip) > 3 holds, implying that 
OV{G,4:,tp) > 0V{G,4:,(p). In order to see that we make a case distinction according 
to the number of neighbours of vertex vi embedded together with vi in the same basis 
subtree. Assume this number is 1 and w.l.o.g. vertex V2 is mapped together with vi to the 
leaves of the same basis subtree, say Ti. Then of course V4, vj and vi are not mapped 
by ip into leaves of Ti. So a{ip) > 3. Moreover, due to the contiguity of ip for at least 
one of the paths {w4,W5,U6}i {^iiViiiVQ}, {^I0 5'?^ii)''^i2} holds that not all of its vertices 
are mapped into the leaves of a common basic subtree. Due to that there is definitely one 
more edge (not incident to vertex vi) whose end-vertices are mapped by ip into leaves 
of different basic subtrees, and hence a[ip) > 4. The other cases where the number of 
neighbours of vi mapped together with vi into the leaves of the same basic subtree is 2 or 
3 can be argued upon analogously^. 



3 A lower bound 

In a DAPT{G,d) with vertex set V{G) of size n, n := |1^(G)|, we have b := d^ leaves, 
where h = [log^^n] is the height of the regular tree. Thus there are i^J^m possible 
arrangements and the complete enumeration becomes inefficient even for very small in- 
stances. Further let us notice that 2m < OV{G, d, (p) < 2hm holds for every arrangement 
(p, where m is the number of edges of the guest graph G. These bounds are due to the 



^In fact we can show that the DAPT is polynomially solvable in the case that the guest graph is an 
extended star as in this example and for some suitable choices of d. In this case the optimal arrangement 
has a particular structure and is in general not contiguous. This and other polynomially solvable special 
cases of the DAPT are discussed in another paper we are working in. 




Figure 3: The guest graph of the DAPT instance in Example^ 




Figure 4: The non- contiguous optimal arrangement which objective value is 28 for the 
DAPT instance in Example\^ 



fact that the distance between any two leaves in a regular tree of height h is between 2 
and 2h. 

Next we introduce the so-called degree lower bound for the DAPT which will be 
also used to evaluate the performance of the heuristics introduced in this paper. We 
adapt an idea used by Petit in [7j for the linear arrangement problem. The idea is the 
construction of locally optimal arrangements for every vertex v of G, i.e. the construction 
of an optimal arrangement of v and its neighbours. Than the contribution of vertex v to 
the objective function value of any feasible solution cannot be larger than the objective 
function value of this locally optimal arrangement divided by 2. 

More precisely, for every v € V{G) we define a new graph G'^ = (y^,E'^) with the 
vertex set Vy := V and the edge set £'(, = {{v,u}: u G r(t>)}. Thus G'^ is a subgraph 
of G containing all vertices of G and just the edges incident to v. Obviously, G'^ is the 
union of a star and some isolated vertices. An optimal arrangement (p^ for DAPT^G'^, d) 
is obtained by placing v on some /ea/, say 61 w.l.o.g. and the other neighbours on the 
leaves 62, • • • , 6i+|r(i;)| one by one, where the canonical order of the leaves is adopted. The 
other vertices of G are arranged arbitrarily on the remaining leaves &2+|r(t))h • • • ) b^h . Let 
OVv denote the objective function value of the above mentioned arrangement for every 
V ^V. It is obvious that DB{G,d) given as below is a lower bound for DAPT{G,d), 
that is 

DB{G, d) = ;^ ^ OK < OV{G, d, <t>) for ah arrangements (j). (3) 

vev 

This bound DG{G, d) is called the degree hound. 

DB{G,d) can be easily computed because OVy can be easily computed, given d and 
the number |r(f)| of neighbours, for all v G y{G). 

Lemma 5 Let G = (V, E) he a star graph with n vertices and 2 < d < n a natural 
number. The optimal value OV of DAPT{G, d) is given as 

( d^ -\\ 

OV = 2 i h n — - — — J , where h = [log^ n] is the height of the host d-regular tree. 

(4) 

Proof. Let v := vi he the central vertex of G with vertex set {vi,V2, . ■ ■ ,Vn}. 
It is clear that the optimal arrangement places the vertices vi,V2, . . . ,Vn into the 
leaves bi,b2,. . . , bn of the d-regular tree of height h, respectively, where the leaves 
are given in the canonical order. Consider a partition of the set of leaves into sets 
Bj = {b is a leaf : dT{bi,b) = 2j} with j = 0, . . . ,h. It is clear that Bq = {bi}, 
Bi = {62, • • • , bd}, CLiT'd hence \Bq\ = 1, \Bi\ = d — 1. 

Generally, for j = 0,1, ... ,h, a d-regular tree of height h contains d -^ d-regular 
subtrees of height j. Clearly one of these subtrees, say Ti contains bi. This subtree has 
in turn d d-regular subtrees of height j — 1 and (only) one of those contains 61 . The set 



Bj consists exactly of the leaves of those d-regular subtrees of height j — 1 of Ti which 
do not contain bi. There are clearly d — 1 such subtrees with d^~^ leaves each. Hence 
\Bj\ = {d- l)d^-^ for all j = l,2,...,h. 

Due to h = [log^n] we have d^^^ < n < d^ and hence the leaves of the basic subtree 
which contains 61 (and thus hosts vi) are all occupied. Consequently the other basic 
subtrees have exactly n — d^~^ > occupied leaves. Thus we get 

h-l h-l 

OV = J2 2j|^jl + 2/i(n - d''-^) = 2{d - 1) ^ jd^-^ + 2h{n - d^-^). 

Using Ylj=i ^"'"^J = "^ 7Tr? ^^ ^^^ ^^^ lemma. D 

By applying Lemma [5] to evaluate OV^ in ^ as the optimal objective function value 
of the DAPT with a guest graph being star graph with \T{v) + 1| vertices we get: 

Theorem 6 Let G = (F, E) be a graph and 2 < d < n a degree of the arrangement tree. 
Then the degree bound is given as 

DB{G,d) = Y, Uv){\T{v)\ + 1) - ^^^ j (5) 

where 

piv):=\logA\r{v)\ + l)]. (6) 

4 Heuristic approaches for the DAPT 

In this section we will introduce some simple greedy heuristics, a construction heuristic 
and two local search heuristics for the DAPT. 

4.1 Simple greedy approaches 

A simple greedy strategy considers the leaves of the guest graph in the canonical order. 
The first leaf is occupied by a vertex selected at random. Then we consider the next leaf 
in the canonical order, place at it "best possible vertex", and repeat this process until 
all vertices of the guest graph have been placed to some leaf. "The best possible vertex" 
means here a vertex which leads to the biggest increase in the objective function value of 
the DAPT. We call this heuristic G2. G2 is a leaf-driven heuristic. Clearly there are also 
vertex-driven greedy algorithms which investigate the vertices in some prespecified order 
and place the current vertex to the "best possible free leaf". Since the vertex-driven 
greedy heuristics we have tested were outperformed by the leaf-driven greedy heuristic 
described above we do not present them in details in this paper. 



The time complexity of G2 is 0(max{(m + ?7-)n,n^ logn}). To see this consider 
first a pre-processing step to compute the distances between all pairs of leaves of the 
arrangement tree in 0(n^ log n) time according to Observation [2j Then n iterations are 
performed to arrange the vertices one at a time. The computation of the increase in the 
objective function value resulting by placing a specific vertex v onto the current leaf takes 
0(|r(t;)|) time per each vertex and hence 0{m) time for all candidate vertices. Selecting 
the best among all candidate vertices takes another 0{n) time. Thus we obtain a time 
complexity of 0{n+m) per iteration which results to 0{{n+m)n) for all iterations and to 
an overall time complexity of 0(max{(7Ti + n)n,n^ log n}) (including the pre-processing 
step). 

We have also tested two very simple search heuristics BFSG and DFSG which order 
the vertices of the guest graph according to breadth-first search or depth-first search, 
respectively, after starting at some prespecified vertex. Then the vertices are places onto 
the leaves in the canonical order, i.e. the i-the vertex according the resulting ordering is 
placed at the i-th leaf, i = 1, 2, . . . , n. 

Of course there are a number of variants of this algorithm. We distinguish different 
implementations for connected and non-connected graphs. In the case of a connected 
guest graph G there is a flexibility in choosing the starting vertex for search algorithm in 
G. Depending on the graph structure the vertex with the highest degree can be chosen. 
Or the algorithm is run for each vertex as starting vertex and then the best obtained 
solution is chosen. 

In the case of non-connected graphs we have to fix the order of the connected compo- 
nents before running the search algorithm for each of them. This can be done in many 
ways, e.g. by considering the connected components in decreasing order of magnitude. 

Clearly, the worst-case time complexity depends on the particular implementation in 
each case. In the case of connected graphs we obtain an 0{rv') algorithm, if the "best" 
starting vertex among all is chosen. In the case of non-connected graphs we obtain the 
same time complexity, if we choose the best starting vertex in each component by running 
the algorithm as many times as the nuinber of vertices for each component. 

4.2 A construction heuristic 

Let us now consider the objective function of the problem from another point of view. Let 
tti, 1 <i <h,he the number of edges of the guest graph G whose endpoints are mapped 
into leaves of T at a distance 2i in the host graph. 
We can state obviously 

OV{G, d, </.) = 2hah + 2{h - l)ah-i + • • • + 2ai , (7) 

where ah + o/i-i + • • • + ai = m and m is the number of edges of the guest graph G. 

Since our aim is to minimise the objective value OV{G, d, (p), we try first to minimise 
the coefficient Uh by partitioning the vertex set V in at most d subsets Vi, 1 < i < d, 

I Dl 

with < |Vi[ < ^. Then each Vi, 1 < i < d, is embedded into the leaves of the 
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corresponding basic subtree, which means that the inequahties {i — l)d-\- 1 < 4>{v) <id 
hold for any v & Vi, 1 < i < d. Among all arrangements of this kind we choose 
one which minimises ah = \{{u,v) € E\u &Vi,v (z Vj,i ^i}|. Then the subproblems 
DAPT{G[Vi],d), 1 < i < d, (where G[Vi] is the subgraph of G induced by the set of 
vertices Vi) are solved in order to determine an arrangement of Vi, 1 < i < d, into the 
leaves of the corresponding basic subtree. 

The problem of partitioning V as described above is strongly related to the so called 
minimum cut problem with hounded set size (MCBSSP) described in next subsection. 
In Subsection 14.2.21 we present an approach to solve the DAPT{G, d) by using the idea 
described above and a heuristic for MCBSSP. 

4.2.1 A related problem (MCBSSP) and some heuristic approaches 

The Minimum Cut Problem with Bounded Set Size (MCBSSP) 

Input: A graph G = iV^E) with n = \V\ and two integers /, u with < / < 

u < n. 

Output: A set X C F with I < \X\ < u such that the cut 

5{X) := {{u,v) G E\u & X,v ^ X} has minimum cardinality. 

MCBSSP is equivalent to the so-called {k,n — k) cut problem (k-{n — k)CP), inves- 
tigated by Feige, Krauthgamer and Nissim [Sj. 

The {k,n — k) cut problem (k-(n — k)CP) 

Input: A graph G = {V, E) with n = \V\ and an integer k, with k < n. 

Output: A partition of F in X, y with \X\ = k, \Y\ = n — k such that the 

cut 

5{X) := {{u,v) E E\u G X, u G Y} has minimum cardinality. 

Indeed the equivalence between MCBSSP and k-{n — k)CP is trivial: an optimal solution 
of MCBSSP in a graph G with input parameters l,u can be obtained by solving 0{n) 
instances of k-{n — k)CP in the same graph G with input parameter k = u,u + 1, . . . ,1. 
On the other hand k-{n — k)CP is just a special case of MCBSSP, when u = I holds. 
k-{n — k)CP is NP-hard for general k as mentioned in Feige et al. [3], a special case of 
it is the minimum bisection problem, see Carey and Johnson [4J. Thus MCBSSP is also 
NP-hard for general / and u and there is no hope to optimally solve it in polynomial 
time (unless P = NP). 

We have considered two heuristic approaches to solve MCBSSP. These will then be 
applied recursively to obtain a heuristic for the DAPT(G, d) as described above. 

The first approach is based on a polynomial time approximation algorithm for k- 
{n — k)CP with an approximation ratio 0(log n) proposed by Feige, Krauthgamer and 
Nissim [3J. (Their algorithm reaches an even better approximation rate for the cases 
k = O(logn) and k = Vtilogn)). So in order to obtain a solution of MCBSSP in the 
graph G with parameters / and u we apply the approach of Feige et al [3| to k-{n — k)CP 

10 



in G with parameter k varying between / and u and then choose a minimum cut among 
the I — u + \ obtained solutions of k-{n — A;)CP. Since u — l<nwe get a polynomial time 
approach for MCBSSP. 

Our second approach for MCBSSP makes use of a simple local search idea. Assume 
that / = u. We randomly partition y in X and V\X, where ^ C X CV and \X\ = I = u. 
We try to decrease the cardinality of the cut |(^(^)| by the following pair-exchange 
approach. Consider an other cut 6{{X \ {u}) U {v}) for each pair {u,v), where u & X 
and V (^ X. Replace X by (X \ {u}) U {v} if 6{{X \ {u}) U {v}) < 6{X) and repeat this 
step until no further improvement of the cardinality of the cut is possible. Then apply 
the above approach to determine a cut 6{X^^') with [X'^^'l = k for any I < k < u and 
choose the best among the cuts 6{X^^'), I < k < u. 

4.2.2 A heuristic for DAPT(G,d) 

Having described the heuristics for MCBSSP let us turn back to the DAPT{G,d). The 
approach is presented in the form of a pseudo code in Algorithm 14.11 and involves the 
heuristic solution of the MCBSSP as a subroutine (see pseudocode line 11). 

We first consider the question of determining the "unused leaves" , i.e. leaves of the 
arrangement tree, into which no nodes of the guest graph are arranged. Based on our 
observations in the context of numerical tests we try to use as few basic subtrees as 
possible to arrange all nodes of the guest graph. Thus we collect the unused b — n 
leaves (recall that b := \B\ is the number of leaves of the host d-regular tree) into as few 
basic subtrees as possible. By considering that each basic subtree has 6i := ^ we mark 



the first L, 



b~n 



b—n 
bi 



basic subtrees as unused 



bi leaves, or equivalently the first 
(see pseudocode lines [7] - [9|) . Then we separate the vertices X which will be placed on 
the leaves &;„„+!, ... 6, , & , i.e. on the leaves of the first used basic subtree, by solving 

uu~r ^ 

MCBSSP with the parameters I := bi — {b — n) mod 6i and u := 6i (see pseudocode 
line llip . This can be done by applying one of the heuristics described in Subsection l4.2.11 



We repeat then this procedure 



6i 



subproblems which are solved 



-1 times to obtain 

recursively (pseudocode line ll2p . The recursion calls will terminate when the height of 
the arrangement tree becomes 1; there an arrangement (f) is selected at random. 



Now let us consider the worst-case time complexity of the described approach. Let 
fc{n) denote the worst-case time complexity of the subroutine which solves MCBSSP 
for a graph with n vertices and any parameters < I < u < n. Since n <b holds for all 
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Input: G = (y, E) undirected graph and positive integer d G N where 2 < d < n; let be 

\V\ = n and T the d-regular arrangement tree with the set of leaves B 
Output: arrangement (/) : V —?■ B 



1: h := [log^n] and b := h'^; 

2 

3 

4 

5 



9 

10: 
11 



if /i = 1 then 

make the arrangement (p at random; 
else 

for i := 1 to d do 
if luu > 2 *hen 

4>'^^{l) := unused, {i — 1)^ <l<i 
h ■=h - ^■ 



d' 



else 

find a minimum cardinality cut X C V{G) in graph G subject to ^ — buu < 
1^1 ^ ^ by solving MCBSSP with parameters I := ^ — luu and u := ^; 
12: solve the problem for the graph G[X] and a d-regular arrangement tree Tx 

which height is h — 1 recursively; let (j)x be the solution of this recursive 
problem; 
compute the inverse function of cpx which we denote 0^ ; 
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14: 
15 
16 

17: 
18 
19 
20 
21 
22 



for j := 1 to ^ do 

rH{^-l)^,+j)■■=rxHj)■, 

end for 

G := G[G\X]; 

I'UU ■ — I'UU (,^ 1^ Ijj 

end if 
end for 

compute the function (j) from the function (j)^^; 
end if 

Algorithm 4.1: Construction heuristic. 



instances, the worst-case time complexity of the whole algorithm is 
4/c(|<<)+/c(|(^-l))+^^+/c(f2)) + 



(8) 



where the lines correspond to the recursion depth. Summarising we get the following 
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worst case time complexity 

h-2 d-2 



Y.d'Y.fc[:^^d-j)]. (9) 

j=0 jr=0 



For some particular heuristic to solve the MCBSSP we can substitute fcin) by a 
precise expression in ([9]). Consider the case of the local search based heuristic described in 
Subsection l4.2.1l When computing the cuts at the first recursion level u — l<2 obviously 
holds. If X is the set of vertices generating the cut, then \X\ < ^ holds. When computing 

the fc-th cut at the first recursion level we have at most 2 d vertex pairs which could 
be exchanged and the cardinality of the cut after the pair-exchange can be computed 
in 0(^ + (d — A:)^) = 0(^(d — A: + 1)) time. So we get a worst-case time complexity of 

O (^ (^^) (^ + ^)) =0(l,(}^^)W-k + l))= 0{{\f {d-k)id-k+l) 
for the A:-th cut in the first level (where the first factor in the above expression accounts 
for the number of k-{n — k)CP to be solved which is at most u — I < 2)- Summarising 
for all cuts of the first level we get 




Now let us consider the recursion. After building the first d — 1 cuts we get d 

h 

d 



subproblems each of them having most -^ vertices. Thus for the whole algorithm we get 



a time complexity K with 

K:=oi^^+d^+d^^^ + --- + d^-^^^^+r!\. (11) 

\dd d d I 

Note that if the height of the arrangement tree is 1, the arrangement (j) can be made 
at random and thus the recursion depth is only h — 2. Using d ^"^ J = b-^^J— — - 



4 

b 

T 



d' 



^h i \d ^ J _ _\_dj_ — jj^2 ^^^ considering b < nd we get 



^- = o{^t0-M^-'^^nYo(u',). (12) 



Now, we can state the following theorem. 



Theorem 7 The Algorithm ^4jJ\ can be implemented with a worst case time complexity 
of 0{ri^d), if the local search approach of Subsection \4^2A is applied to solve MCBSSP. 
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In fact the quality of this construction heuristic depends significantly on the quality 
of the heuristic used to solve MCBSSP. However, even if we were able to solve MCB- 
SSP to optimality, the construction heuristic would not necessarily compute an optimal 
arrangement. As an example consider DAPT{G, 2) with guest graph G as shown in 
Figure m Figure [6] shows an arrangement obtained by the construction heuristic, where 
MCBSSP was always solved to optimality during the algorithm. This arrangement is 
not optimal; a strictly better arrangement is shown in Figure [7] (this is actually an opti- 
mal arrangement). The reason for this behaviour relies on the fact that minimising the 
coefficients a^, i = 1, 2, . . . , /i, starting with a^ and proceeding in the above order, does 
not necessarily lead to a minimum value of OV{G^ d, </>), see ([7j). 

In our computational experiment we observed that the construction heuristic which 
involved the pair-exchange approach to solve MCBSSP outperforms the heuristic which 
involves the approach of Feige et al. [3]. Therefore in Section [6] we just report on the 
performance of the more successful algorithm denoted by CHLS, see also Section [6741 




Figure 5: Graph. 




Figure 6: A non-optimal arrangement 4> with OV{G,2,(f)) = 26 for G in Figurel^ 

4.3 Local search approaches 

In this paragraph we propose two different local search heuristics for the DAPT. They 
can be used separately or also combined as described below. 
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Figure 7: An optimal arrangement cp with 0V{G,2,(p) = 24 for G in Figure\^ 

4.3.1 The pair-exchange heuristic 

The algorithm starts with an arbitrary arrangement (p (it can be a random arrangement 
or an arrangement obtained by applying some other heuristic) and tries to improve 
the objective function value by performing so-called pair-exchanges. More precisely the 
algorithm fixes an ordering of the pairs of vertices {vi,Vj) G V{G) x V{G) with Vi ^ Vj 
and checks whether a pair {vi,Vj) exists such that OV{G,d,(j)') < OV{G,d,(l)), where (j)' 
is obtained from (j) by applying a pair-exchange: 

{(j){vj) a k = i 
(Pivi) iik = j forA;E{l,...,n}. (13) 

(pivk) iik^{ij} 

If such a pair {vi,Vj) of vertices whose exchange improves the objective function value can 
be found, then cp is substituted by cp' and the procedure is iteratively repeated. Otherwise 
the algorithm terminates and outputs the current arrangement. Note that this approach 
would keep unchanged the set of unused leaves. In order to be able to vary it we work with 
an extended guest graph G' = (V' , E') with vertex set V = V U {vn+i-, ■ ■ ■ vi,} and edge 
set E' = E, where V is the vertex set of the original guest graph G. The new guest graph 
has as many vertices as the number of leaves of the host graph. Since the vertices Wn+i, 
. . ., vi) oiV'\V are isolated vertices, OV{G, d, <p) = OV{G' , d, cp') obviously holds for all 
arrangements cp : V —i' B and all arrangements cp' : V' ^f B such that (p{v) = (p'{v) for 
all V ^V. Thus we can solve DAPT{G\d) instead of solving DAPT{G,d); an optimal 
solution 0* of DAPT{G,d) is obtained from an optimal solution cp'^, of DAPT{G' ,d) 
by setting (p^{v) = (p'^:{v), \/v € V. Note that, however, if the starting arrangement is 
contiguous, then applying the pair-exchange to DAPT{G' ,d) instead of DAPT{G,d) 
can not generate any variation in the set of used leaves. The reason is that 

dT{bi,bj)<dT{bi,bk) (14) 

holds for al\l<i<j<k<b and thus a pair-exchange which arranges an isolated 
vertex v' € V'\V between some pair of eventually connected vertices can never improve 
the objective function value. 
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Theorem 8 The pair-exchange heuristic for the DAPT{G, d) can he implemented with 
time complexity 0(n^(i^mmin{?7z,n}(logn)), where n is the number of vertices and m is 
the number of edges in G. If the starting arrangement is contiguous, then the heuristic 
can be implemented in 0(n^mmin(m,n)(logn)) time. 

Proof. There are 0(6^) pairs of vertices in the graph G' . Since 2m < OV{G,d,(p) < 
2hm holds for every arrangement (p, we can make at most 0{2hm — 2m) = 0{hm) = 
0{{).ogh)m) = 0{\og{nd)m) = O {{log n + log d)m) = 0{mlogn) improvements of the 
objective function value (if d is considered to he a constant and by using b < nd). 

Consider that the pairwise distances between all pairs of leaves in the arrangement 
tree can he computed in 0(6^ log n) = 0(n^d^ log n) time in a pre-processing step, see 
Observation \^ In order to update the objective function value of an arrangement after 
a pair- exchange of vertices Vi and Vj which transforms the current arrangement (p to 
the arrangement (j)' as in il3\) . the length of the path between (p{vi) (4>{vj)) and <j){v) 
is substituted by the length of the corresponding path between (j)'{vi) (4>'{vj)) and (j)'{v), 
for all neighbours v of Vi (vj). Since the vertices which exchange position have in total 
0(inin{7Ti, n}) neighbours, the objective function after a (candidate) pair-exchange can 
be updated in 0(min{m-,n}) time. With at most 0{b'^) (candidate) pair- exchanges to he 
performed in each iteration and at most 0{mlogn) iterations, the overall time complexity 
of the algorithm amounts to 0(6^ min{r?2, njmlogn) = 0(n^(i^rn-min{m,n}logn). 

// the starting arrangement is contiguous, then just 0(n?) pairwise distances need 
to be computed in the pre-processing step and the overall time complexity amounts to 
0(n^7nmin{m, njlogn). D 

Clearly, we can also fix an ordering of the pairs of leaves and exchange the vertices 
arranged at some pair of leaves (if any), in this ordering. One would obtain a similar time 
complexity as in the general case of Theorem [HI We refer to these heuristics as vertex- 
based pair-exchange heuristic and leaf-based pair-exchange heuristic, respectively. Our 
computational experiments have shown that the vertex-based pair-exchange heuristic 
generally outperforms the leaf-based pair-exchange heuristic. For this reason we only 
report about the performance of the vertex-based pair-exchange heuristic (abbreviated 
by PEHVNA) in Section [3 

4.3.2 The shift-flip heuristic 

The last heuristic we discuss is the shift-hip heuristic. First, we need two definitions. 

Definition 9 {Fhp) Let G = {V,E) be an undirected guest graph with \V\ = n, T a 
d-regular tree, with 2 < d < n, and let B be the set of leaves ofT. Let cj) : V ^>- B he an 
arrangement. Further, let e, 5, /, r € N U {0}, he parameters with < e < h, 1 < g < d^, 
1 < I < r < d. Finally let f be a bijection f : B ^- B defined as follows: 

l)d^-(e+i) +ti,l<ti< d'^-(^+i) 
fib^) = { b^,.,a.iu^-i.+^)+u for i = Ho) + {r - l)d^-^^+^) +ti,l<ti< d^-(^+i) , 
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^A(g)+(r-l)dft-('=+i)+ti 


for i = A(5) + (/ 


^A{s)+(/-l)d'»-(=+i)+t, 


for i = A{g) + (r 


hi 


otherwise 



(15) 



where A{g) := {g — l)d . The arrangement (pj : V ^ B where (f)f = f o (p is a flip of 
tlie arrangement (p. We say that we Eip the arrangement (j) at the l-th and r-th d-regular 
subtrees of the g-th node in level e. 

In a more descriptive explanation a flip consists of intercfianging the preimages of 
the leaves of two d-regular subtrees of the arrangement tree which have the same height 
and whose roots have a common father vertex, while preserving the order of the leaves 
in each of the two interchanged subtrees. More precisely we consider the vertices of the 
d-regular tree as being partitioned into levels, the root having level 0, its d sons having 
level 1 and so on, to end up with the leaves at level h — 1. In Definition [9] we consider the 
g-th vertex in level e and the indices / and r of two sons of that vertex. The successors 
of each of those suns build a d-regular subtree of height /i — (e + 1), respectively. The flip 
operation interchanges exactly the preimages of the leaves of these two d-regular subtrees 
by preserving in each subtree the order of the leaves induced by the canonical order of 
the leaves in T. 

For an illustration consider an instance DAPT{G, d) with guest graph G given in 
Figure M and d = 3. 



Figure 8: A guest graph. 

Consider further an arrangement represented in Figure [9l each filled leaf contains the 
index of the vertex of G mapped into that leaf. 




Figure 9: An arrangement (j) with OV{G, 3, 



32 for G in Figure 



In Figure [10] we see the flip obtained from the arrangement represented in Figure [9] 
with parameters e = l, g = 2, 1 = 2 and r = 3. Note that flipping does not change the 
objective function value of the arrangement. 
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Figure 10: A flip of the arrangement shown in Figure 0, for the guest graph shown in 
FigurelM and d = 3. The parameters of the flip are e = l, g = 2,l = 2 and r = 3. The 
objective value of the flipped arrangement remains unchanged and equals 32. 

Proposition 10 Let G = {V, E) he an undirected guest graph with \V\ = n,T a d-regular 
tree, with 2 < d < n and let B he the set of leaves of T. Further, let e, 5, /, r € N U {0}, 
he parameters with < e < h, 1 < g < d'^, 1 < I < r < d. Let f he a hijective 
function f : B —^ B defined as in Definition For any arrangement (p : V ^ B 
and the corresponding flip (j)f : V —^ B of the arrangement (j), (j)f = f o (j), the equality 
OV{G,d,^f) = OV{G,d,^) holds. 

Proof. To prove the statement we make use of Observation [2j For Vi €z V, i = 1,2, . . . ,n, 
let us denote by p{i), Pf{i) the indices of the leaves <j){vi), (t)f{yi) of T in the canonical 
ordering, respectively. We clearly have p{i),Pf{i) £ 1, 2, ... , d'^, for all i = 1, 2, . . . , n. 
According to Observation [2] we get the following expressions for the objective function 
values of (f> and (pf: 



OViG,d,. 



^ dT{4>{Vi),(t){Vj)) 



(16) 



{vi,vj}eE{G) 



{vuVj}£EiG) 



yj 2argmin <k£{l,2,...,h}: 



p{i 



d^ 



PJj) - 1 
d^ 



and 

OV{G,d,cPf) 



^ dT{(t}f{Vi),(t)f{vj)) 



{v„Vj}&E{G) 



{v„v,}&E(G) 
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y^ 2argmin<^ A; G {1,2, ... ,/i}: 



Pf{i) 



d^ 



Pfij) - 1 
d^ 



Consider the index p of an arbitrary leaf hp of T (in the canonical order) written as 
p = {u — l)d^~'^ + (s — l)d^~^^~^^' + 1 for some natural numbers 1 < u < d^, 1 < s < d 
and 1 < t < d ''^"^^^ u represents the index of the unique node x at level e which is an 
ancestor of bp, s represents the index of the d-regular subtree Ti of height /i — (e + 1) 
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hanging on x and t represents the index of bp in Ti according to the canonical order of 
the leaves of Ti induced by the canonical order of the leaves of T. Then the following 
equality holds 



p — 1 



{u - l)d^-^-'' + (s - l)d''-(e+i)-fc + ^tz-i J if fc < /i - (e + 1) 
(u — 1) \i k = h — e 

a k > h — e 
iik = h-(e + l) 



, (18) 



u-l)d + {s-l) 



for any 1 < u < d'^, any 1 < s < d and any 1 < t < d^ (^=+1). Notice that according 
to Definition El (/)(t;i) / (j)f{vi) holds, only if p(z) = A{g) + (/ - l)(i^~(^+^) + U or p{i) = 
A{g) + (r — l)d'^~^'^~^^^ + ij with some 1 < ti < d'^~^'^^^\ Moreover, the following two 
implications hold for tj = 1, 2, . . . , d'^~^'^~^^': 



p{i) = A{g) + {l- l)(i''-(^+i) + ti implies 



Pf[ 



A{g) + (r - l)d 



h-{e+l) 



+ ti, (19) 



pii) = A{g) + (r - l)d^-(^+^) + U implies pf{i) = A{g) + {I - l)d^-(^+i) + U. (20) 



Consider now an edge {vi,Vj) with (/>(fi) 7^ </'/(^'i) or (/>(i'j) 7^ (l)f{vj), which is equiv- 
alent to p{i) 7^ Pf{i) or p(j) 7^ Pfij)- There are two cases: (I) p{i) 7^ Pf{i) and 
p(j) 7^ PfiJ), or (II) just one of the inequalities p(i) 7^ Pfii), pij) 7^ P/(j) holds. 

Case I. In this case one of the following cases can happen: 

Case la. p{i) = A{g) + {I - l)d^-(^+^^ + U and p{j) = A{g) + (/ - l)d^~^^+^) + tj, or 

Case lb. p{i) = A{g) + (r - l)(i^-('^+i) + ti and p[j) = A{g) + (r - l)(i'^-(^+i) + tj, or 

Case Ic. p{i) = A{g) + {I - l)d^-^^^^) + ti and p(i) = A{g) + (r - l)d^-(^+i) + tj, or 

Case Id. p(z) = A{g) + (r - l)d'^-(<^+i) + ti and p(j) = A(c/) + (/ - l)d'^-(<=+i) + t^-. 

In Case Ic and in Case Id we get d{(j){i), 4>{j)) = d{(j)f{i), 4>f{j)) = 2{h — e) by applying 
p^ and considering ([19]), ([20]) . In Case la and in Case lb we get 



d(<A(i),0(i))=d(0/(i),0/(j)) = 

2min<^ h - {e + l),argmin|A; e {1,2,/i - (e + 2)}: 



Zi 1 tj 



d^ 



d^ 



'-}]^ 
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Case II. Assume w.l.o.g. that p(i) = {g - l)d^~'' + {I - l)(i^"(^+^^ + U and let 
p{j) = {u — l)d^~^ + (s — l)dP'~^^^^' + tj, where g ^ u oi s ^ {h^}- Clearly Pf{i) = 
{g - l)d^-^ + (r - l)d'^-(^+i) + ti and pf{j) = p{j) = {u - l)d^-^ + {s - l)d'^~^^+^^ + tj. 
If u = g and s {l,r}, then (|18p together with Observation [2] implies dT{(f>{i)4'{j)) = 
dT{(t>f{i)(t>f{j)) = 2{h-e). 



Otherwise, \i u ^ g, then (llSp implies 



p(i) 



d'^ 



/ 



p(i) 



d'^ 



p/( 



d^ 



p(i) 


-1 


= 


pj_ 




y^ 


p(i) - 
d^ 


1 


= 



for all A; > /i — e and 
PfU) - 1 



d'^ 



for all k < h — e. 



which together with Observation [2] implies then dT{(pii),4'iJ)) = dT{4>f{i)-,4>f{3))- 
Thus dT{(t){i)4>{3)) = dTi4'f{i)4'f{J)) for aiiy edge {vi,Vj) £ E. Therefore the right-hand 
sides of the equations (fT6|) and pT|) are equal and OV{G, d, cp) = OG{G, d,(j)f). D 

Definition 11 {shift) Let G = {V^E) he an undirected guest graph with \V\ = n and 
T a d-regular arrangement tree with 2 < d < n, set of leaves B and number of leaves 
b = \B\. Let (j) : V ^ B he an arrangement. Further, let k a N be an integer. An 
arrangement cpk with 



<j)k{v) := 6(((j_i)+fc) mod b)+i , where (piv) = h , 
is a shift of the arrangement (p. We say that we shift the arrangement (p by k. 
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The idea of the shift-flip heuristic is fairly simple. For a given arrangement <p we 
find out a 1 < k < b which minimises the objective function value OV{G,d,(pk)- There 
are two possibilities to define the shift step. In the first variant we apply the shift by 
k defined as above only if it implies an improvement of the objective function value, 
i.e. OV{G,d,(pk) < OV{G,d,(p), and substitute then the current arrangement (p by the 
improved one (pk- In the second variant we also accept an arrangement which keeps the 
objective function value unchanged. If no such an arrangement can be found, then a 
further fiip is performed. Both approaches proceed in the next iteration by applying 
a random fiip to the current arrangement and so on until a termination criterion is 
satisfied. Both variants of the heuristic output the best arrangement found during the 
search. We report about the performance of the second variant because this variant 
seems to outperform the first one. 

Of course there are a number of possibilities to define a terminating criterion. It can 
be a run time bound which defines the maximum length of a time interval the algorithm 
is allowed to run without doing an improvement. Or it can be a bound on the overall 
number of fiip and shift steps performed without improving the objective function value. 

Both variants of the shift-fiip heuristic (SF) can be combined with the pair-exchange 
heuristic (PE). Since the search neighbourhoods of the two heuristics are significantly 
different, it is possible to escape from the local minima of SF by just applying a search 
in the PE neighbourhood and vice-versa. 
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5 Test instances 

We test and compare the above described heuristics on some famihes of test in- 
stances. To the best of our knowledge there are no standard test instances 
for this problem, so we have generated some test instances ourselves. We 
introduce the following families of test instances which are also available at 



http : //www . opt .math . tu-graz . ac . at/~cela/public . htm 



5.1 Test instances solvable by complete enumeration 

The guest graphs of these instances are marked by the prefix "CEJ'. The first graph in 
this category CE_sample corresponds to the graph in Figure [H Further we consider 2 
(thin) graphs CE_thin7 {n = 7, m = 7) and CE_thinlO (n = 10, m = 11) with 7 and 
10 vertices, respectively. We generate test instances with guest graph CE_thin7 and all 
possible values of d, 2 < d < 7. With the guest graph CE_thinlO we generate instances 
with d = 2 and d = 4. Further we consider some analogous instances with denser guest 
graphs: CE_dense7 {n = 7, m = 14) and CE_denselO (n = 10, m = 26) with 7 and 10 
vertices. Finally, we consider a 3 x 3 mesh CE_mesh9 and d = 2, d = 3 and d = A. 

For this family of test instances the precise values of d were chosen so as to be able 
to solve these instances by complete enumeration within a prespecified time limit, see 
Section [6l 

5.2 Test instances with known optimal solution 

These instances are special cases of the DAPT which can be solved by a polynomial 
time algorithm, see [UdO]. The guest graphs of these instances are marked by the prefix 
"SCJ'. Unless the special case involves a particular choice of d, we use d = 2 and d = 7 
for all considered guest graphs. We consider instances of following types: 

• Instances for which d = n — 1 where n is the number of vertices of the guest 
graph. We use 3 guest graphs generated at random for this class of instances: 
SC-randomBS, SC-randomSO and SC-random75. These graphs have the same num- 
ber of vertices, n = 500, and in each of them any pair of non-equal vertices build an 
edge independently at random with probability 0.25, 0.50 and 0.75, respectively. 

• Instances whose guest graphs build a star, that is they consist just of a central ver- 
tex connected by an edge to all other vertices of the graph. The concrete graphs 
are SCstarbO, SC.star500 and 5C_siarlOOO with 50, 500 and 1000 vertices, re- 
spectively. 

• The guest graph in Figure [3] and the choice d = 4. This guest graph is an extended 
star and this instance is referred to as SC-extStar. 
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• Instances whose guest graphs build a d-regular tree. We denote these guest 
graphs/instances by SC-treeDGxHy where x = d holds and y is the height of 
the tree. 

• Instances whose guest graphs build of a path. We denote these guest graphs by 
SC.pathSO, SC.pathSOO and SC.pathlOOO. They have 50, 500 and 1000 vertices, 
respectively. 

• Instances whose guest graphs build a simple cycle. We created 3 graphs of this type: 
SC_simpleCycle50, SCsimpleCycleSOO and SCsimpleCyclelOOO with 50, 500 and 
1000 vertices, respectively. 

5.3 Randomly generated test instances 

The guest graphs of these instances are marked by the prefix "RGJ'. This instances 
are generated in the same way as the instances SC-random25, SC-random50 and 
SC-random75. All guest graphs in this class of instances have 500 vertices and the pairs 
of vertices are present as edges in the graphs randomly and independently with the same 
constant probability, say ■^. For each x two random graphs are constructed as above 
and are denoted by RG-randomAx and RG-randomBx. The degree of the regularity of 
the host tree is set to d = 2 and d = 7. 

5.4 Instances with graphs taken from Petit [8] 

The guest graphs of these instances are marked by the prefix "Pet03_". These graphs 
were used in j8j to test some heuristics for the linear arrangement problem (LAP), a 
problem related to the DAPT as explained in Section [TJ Also in this family of instances 
we use d = 2 and d = 7. This choice of the parameter d is motivated by the goal of 
comparing the behaviour of the proposed heuristics when a smaller and a larger value of 
the parameter d are considered (d = 2 and d = 7). 

6 Numerical results 

The results of all numerical tests are summarised in the tables in Appendix. We group 
the test instances described in Section [5] in three groups: instances solvable by complete 
enumeration, polomially solvable instances and the rest. Table [2] reports on instances 
which could be solved to optimality by complete enumeration on the following computer 
in 1 week: HP Compaq nx7400, 32 bit Intel processor (Intel® Centrino® Duo T2250 1.73 
GHz), running in Ubuntu (Linux). Tab el [3] summarises the results for the instances which 
are solvable to optimality in polynomial time. Table H] summarises the computational 
results obtained for the remaining instances. 
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In order to compare the quality of the proposed heuristics we define a quality quotient 
as fohows 

where J^ denotes a set of test instances DAPT{G, d) with guest graph G and degree 
of regularity d. J^ denotes a set of heuristics, HE{G, d) stays for the objective value 
obtained from the heuristic HE for the instance DAPT{G,d), DG{G,d) represents the 
degree bound for this instance and OS{G, d) stays for the objective function value of an 
optimal solution. We set OS{G, d) = if the objective value of an optimal solution is 
unknown. We also write q{G,d,Jif) for q{J^,J^) \iJ = {DAPT{G,d)}. 

We evaluate also the so-called success factor which for a certain group of instances 
and a certain heuristic gives the proportion of instances for which the considered heuristic 
computes the best known solution. 

6.1 Results on test instances solvable by complete enumeration 

Let us first consider Table [2j All instances of this class are very small (in fact, they have 
only 5-10 vertices), and thus most of the heuristics were able to return an optimal 
solution for many instances. The success factors for the instances of this group are 
summarised in Figure [11] (the acronyms are listed on the last page) . The DB entry shows 
us the proportion of the test instances whose optimal objective function value equals the 
degree bound. The degree bound coincides with the optimal objective function value 
only in the special case d = n. Notice that in this case all arrangements yield the same 
objective function value and the DAPT is trivial. 
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Figure 11: Success factors for the instances solved by complete enumeration. 



6.2 Results on test instances solvable in polynomial time 

Table [3] is related to the instances which can be solved by a polynomial time algorithm. 
This group of instances is divided into four parts as follows. 
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Figure 12: Success factors for the instances with d = n — 1. 
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Figure 13: Success factors for the instance with d = 4 and the guest graph of Figure\^ 
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Figure 14: Success factors for the instances with a d-regular tree as a guest graph. 
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(i) Instances for which the equahty d = n — 1 holds. The corresponding success factors 
are given in Figure [121 It is interesting that CHLS yields the optimal solution 
for any instance of this group, which of course does not hold for all such DAPT 
instances in general, cf. e.g. Figure [HI 

(ii) Instances whose guest graph is a star, a simple path or a simple cycle. Most of the 
heuristics return an optimal solution. 

(iii) The instance with the guest graph of Figure [3] and d = A. The corresponding 
success factors are given in Figure [131 Note that neither the lower bound nor any 
heuristics is able to reach the optimum. Note also that some heuristics can generate 
a non-continuous arrangement. In our implementation they are RAM, CHLS and 
SFHWI. 

(iv) Instances with a d-regular tree as a guest graph. No heuristic is able to return an 
optimal arrangement for these instances and TFSG performs mostly better than 
CHLS. The quality of the solutions is q{^ , J^) ~ 1.18. The corresponding success 
factors are given in Figure [TH 

6.3 Results on test instances with unknown optimal solution 

Let us now consider the test instances with unknown optimal solution, i.e. the optimal 
solution of this instances is not obtained by complete enumeration and it is not known 
whether it can be computed in polynomial time, see Table [H The corresponding success 
factors are given in Figure [T5l 
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Figure 15: Success factors for the instances with an unknown optimum. 

In the following we make some remarks on the particular classes of instances from 
this group. 

Consider first the randomly generated instances (with prefix RG-). For all these 
instances CHLS outperforms the other heuristics. We also observe that for any fixed 
d the quotients q{J^, J^) are better for denser graphs. The overall quality quotient for 
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all these instances is g(j^,^) « 1.21. The quality quotient is better if d = 2; we get 
g(j^,^) ?» 1.20 over the instances J with d = 2 and q{J ,M') « 1.23 for the other 
instances of this group. 

For the random instances {RG-randomAx and RG-randomBx) we also observe 
an improvement of the quality quotient depending on the increasing expected den- 
sity of the guest graph. Figure [TU] shows the values of quality quotient computed 
for each pair of instances with guest graphs RG^randomAx and RG-randomBx, for 
X G {5,15,25,25,45,55,65,75,85,95}, and d = 2 or d = 7 respectively. Clearly x 
represents the expected density of graphs generated as described in Section | 



Next consider the instances with guest graphs taken from Petit |8]. Let us notice 
that we have not considered the guest graphs PetOS-crack with d = 2 and have also 
excluded Pet03-wave and Pet03_small as guest graphs from our tests. The reason is the 
big size of the guest graphs for the first two cases and the obtained solution by complete 
enumeration in the third case. The quality quotient is q{,^,^) ~ 1.84 for this group 
of instances. For d = 2 we get q{J^,^) !^ 1.93 and for d = 7 we get q{J^,J^) ^ 1.76. 
Note that the quality quotient is worse for these instances than for the RG- instances. 

A special behaviour could be observed on following test instances: 

• The guest graph is given by Pet03-hcl0 and d = 2. The underlying graph corre- 
sponds to a 10-hypercube. Five heuristics yield solutions with the same objective 
function value which is the best know so far. It is worth of investigating whether 
this objective function value is optimal. 

• The guest graph is given by PetOS-bintreelO (a binary tree of height 10) and d = 7. 
This problem is polynomially solvable in the case that d = 2 [Ij. For d ^ 2 
the computational complexity of this problem is still open. We observe that TFSG 
performs better than GHLS for both instances with the guest graph PetOS-hintreelO 
and d = 7 or d = 2, respectively. 



6.4 Performance of the construction heuristic 

In Tables El [3] and H] only the variant of the construction heuristic which uses the sim- 
ple local search idea (see Section I4.2.ip to solve MCBSSP is included. This strategy 
outperforms the other one which uses the algorithm proposed by Feige, Krauthgamer 
and Nissim [3J as a subroutine to solve MCBSSP. Table [D provides some results on the 
comparison of the construction heuristic involving both approaches to solve MCBSSP, 
respectively. In this table there is only one instance for which the involvement of the al- 
gorithm of Feige et al. yields better results. The guest graph of this instance is 2-regular 
tree with d = 2, hence this is an instance of a special case of the DAPT solvable in 
polynomial time, see [I]. 
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Figure 16: Progress of the quality quotient q{cy,J^) as a function of the expected density 
of the randomly generated guest graphs. 



Table [T] - Two strategies in our construction heuristic 


Graph 


d 


n 


m 


OS 


LS 


FKN 


SC_random50 

SC_treeDG2H8 

SC_treeDG3H6 


499 
2 

3 


500 

511 

1093 


62468 

510 

1092 


125374 

2434 
3926 


125374 

3466 
5042 


125374 
3188 

5234 


RG_randoniA5 

RG_randoniA55 

RG_randomA95 


2 
2 

2 


500 
500 
500 


6126 

68320 

118499 


- 


86628 
1074734 
1893354 


94660 
1091954 
1898482 


RG_randomA5 

RG_randomA55 

RG_randomA95 


7 
7 

7 


500 
500 
500 


6126 

68320 

118499 


- 


36392 
446038 

785158 


39344 
452220 
787360 


foil 


ow-up on the next page . . . 
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Tabled]- 


- Two strateg 


ies in our construction heuristic 




Graph 


d 


n 


m 


OS 


LS 


FKN 


Pet03_randomAl 


2 


1000 


4974 


- 


71874 


81264 


Pet03_hcl0 


2 


1024 


5120 


- 


56320 


56320 


Pet03_cly 


2 


828 


1749 


- 


16884 


21454 


Pet03_gd95c 


2 


62 


144 


- 


866 


1044 


Pet03_randomAl 


7 


1000 


4974 


- 


29574 


33980 


Pet03_hcl0 


7 


1024 


5120 


- 


25892 


27580 


Pet03_cly 


7 


828 


1749 


- 


7508 


9016 


Pet03_gd95c 


7 


62 


144 


- 


410 


500 



Table 1: A comparison of the two approaches used to solve 
MCBSSP as a subroutine in the construction heuristic: the 
local search idea (LS) and the algorithm proposed by Feige et 
al. ig/ (FKN), see Section\4^ 



7 Conclusions and outlook 

In this paper we deal with the data arrangement problem on regular trees DAPT, identify 
some basic properties and introduce heuristic approaches for this problem. We provide 
a comparative analysis of the proposed heuristics based on a set of test instances we 
have generated. To the best of our knowledge no sources of literature dealing with 
heuristic approaches for the DAPT are available. So there is no possibility to test the 
performance of the proposed heuristics on already known benchmark instances and nei- 
ther to compare the proposed heuristics to already existing approaches in the literature. 
However we make use of test instances available in Petit [8] for a related problem, the 
linear arrangement problem, and use these graphs as a guest graph in our test instances. 
We have summarised the generated test instances in a library which is available at 
[http : //www ■ opt ■ math . tu-graz . ac . at /~c el a/public . htm , 

There is plenty of room for further research on this topic in the future. Most of 
the heuristics we propose are basis approaches which can be well combined with one an- 
other. Especially we expect a significant performance improvement if the two local search 
heuristics we propose are combined in order to may escape form the local minima of our 
neighbourhood by making a jump in the other neighbourhood. Also in the construction 
heuristic there is room for improvement, especially as far as the subroutine used to solve 
MCBSSP is concerned. Since this problem has been investigated to some extent in the 
literature there is hope for appropriate approaches to make use of in the construction 
heuristic. Another aspect which could be considered is an alternative handling of the 
unused leaves. 
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Table [2]- summary for the instances solved by the complete 


enumeration 






Graph 


d 


n 


m 


OS 


DB 


NAM 


RAM 


RCAM 


G2 


BFSG 


TFSG 


CHLS 


PEHVNA 


SFHWI 


CE_sample 


3 


5 


7 


20 


18 


22 


20 


20 


22 


20 


20 


20 


20 


20 


CE.thin? 


2 


7 


2 


24 


21 


26 


24 


24 


26 


24 


24 


24 


26 


24 


CE_thin7 


3 


7 


2 


18 


16 


20 


18 


18 


20 


18 


18 


18 


18 


18 


CE.thin? 


4 


7 


2 


16 


14 


16 


18 


16 


16 


16 


16 


16 


16 


16 


CE.thin? 


5 


7 


2 


16 


14 


18 


18 


16 


18 


18 


18 


16 


18 


16 


CE.thin? 


6 


7 


2 


16 


14 


18 


18 


16 


18 


16 


16 


16 


16 


16 


CE.thin? 


? 


7 


2 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


14 


CE.densc7 


2 


7 


14 


62 


56 


64 


62 


62 


64 


62 


66 


62 


62 


62 


CE.dense7 


3 


7 


14 


44 


42 


44 


44 


44 


44 


44 


46 


44 


44 


44 


CE.dense? 


4 


7 


14 


40 


35 


42 


40 


40 


42 


40 


44 


40 


40 


40 


CE.dense? 


5 


7 


14 


38 


30 


40 


40 


38 


38 


38 


38 


38 


38 


38 


CE.dense7 


6 


7 


14 


34 


28 


36 


38 


34 


34 


34 


34 


34 


34 


34 


CE.dense? 


? 


7 


14 


28 


28 


28 


28 


28 


28 


28 


28 


28 


28 


28 


CEjneshQ 


2 


9 


12 


54 


40 


58 


62 


56 


58 


62 


54 


54 


54 


54 


CEjneshg 


3 


9 


12 


36 


30 


36 


36 


36 


36 


36 


36 


36 


36 


36 


CEjneshQ 


4 


9 


12 


34 


25 


36 


36 


34 


36 


36 


34 


34 


34 


34 


CE.thinlO 


2 


10 


11 


46 


34 


62 


54 


48 


60 


56 


46 


46 


48 


48 


CE.thinlO 


4 


10 


11 


30 


22 


34 


32 


30 


34 


32 


30 


30 


30 


30 


CE.denselO 


2 


10 


26 


134 


118 


140 


150 


134 


136 


134 


140 


134 


134 


136 


CE.denselO 


4 


10 


26 


80 


74 


82 


86 


80 


80 


80 


86 


80 


82 


80 


PetOS.small 


2 


5 


8 


34 


28 


36 


34 


34 


34 


34 


34 


34 


34 


34 


PctOa.small 


3 


5 


8 


24 


22 


26 


24 


24 


24 


24 


24 


24 


24 


24 



Table 2: Summary for the instances solved by the complete enumeration. 



Tabic [3]- summary for the instances solved by a polynomial time algorithm | 


Graph 


d 


n 


m 


OS 


DB 


NAM 


RAM 


RCAM 


G2 


BFSG 


TFSG 


CHLS 


PEHVNA 


SFHWI 


SC.random25 


499 


500 


31239 


62644 


62478 


62720 


124774 


62740 


62644 


62644 


62704 


62644 


62706 


62706 


SC.random50 


499 


500 


62468 


125374 


124936 


125438 


249544 


125410 


125374 


125380 


125418 


125374 


125420 


125438 


SC.random?5 


499 


500 


93548 


187784 


187096 


187832 


373724 


187822 


187784 


187796 


187840 


187784 


187824 


187832 


SC.star50 


2 


50 


49 


474 


286 


474 


478 


474 


474 


474 


474 


474 


474 


474 


SC.star500 


2 


500 


499 


7978 


4488 


7978 


7980 


7978 


7978 


7978 


7978 


7978 


7978 


7978 


SCstarlOOO 


2 


1000 


999 


17954 


9976 


17954 


17968 


17954 


17954 


17954 


17954 


17954 


17954 


17954 
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Graph 


d 


n 


m 


OS 


DB 


NAM 


RAM 


RCAM 


G2 


BFSG 


TFSG 


CHLS 


PEHVNA 


SFHWI 


SC_star50 


7 


SO 


49 


186 


142 


186 


262 


186 


186 


186 


186 


186 


186 


186 


SC-StarSOO 


7 


SOD 


499 


3200 


2099 


3200 


3770 


3200 


3200 


3200 


3200 


3200 


3200 


3200 


SCstarlOOO 


7 


1000 


999 


7200 


4599 


7200 


7600 


7200 


7200 


7200 


7200 


7200 


7200 


7200 


SC.extStar 


4 


12 


11 


28 


23 


30 


32 


30 


32 


32 


30 


30 


30 


30 


SC.treeDG2H8 


2 


Sll 


SIO 


2434 


1529 


8176 


7968 


7982 


8176 


7680 


2746 


3466 


4878 


6426 


PctOS.bintreclO 


2 


1023 


1022 


4904 


3065 


18414 


18146 


18118 


18414 


17410 


5618 


7072 


10696 


1S030 


SC.trccDG2H10 


2 


2047 


2046 


9850 


6137 


40940 


45S08 


40S42 


40940 


38914 


11418 


16026 


23S66 


34938 


SC.trecDG2Hll 


2 


409S 


4094 


19744 


12281 


90090 


89S68 


89S28 


90090 


86020 


23154 


33748 


50826 


80064 


SC.treeDG2H12 


2 


8191 


8190 


39S38 


24S69 


196584 


19S714 


195668 


196584 


188420 


46930 


71666 


109504 


174972 


SC.treeDG3H5 


3 


364 


363 


1296 


967 


3640 


3870 


3642 


3640 


3640 


1524 


1472 


2376 


2778 


SC.treeDG3H6 


3 


1093 


1092 


3926 


2911 


13116 


14040 


13214 


13116 


13116 


4772 


5042 


8380 


11314 


SC.treeDG4H4 


4 


341 


340 


1058 


849 


2728 


3102 


2738 


2728 


2728 


1288 


1328 


1584 


2144 


SC.treeDG4H5 


4 


136S 


1364 


4272 


3409 


1365 


15310 


13884 


136S0 


13650 


5320 


5412 


7626 


11820 


SC.treeDG8H3 


8 


S8S 


S84 


1472 


1313 


3510 


4426 


3S26 


3510 


3510 


1956 


1808 


1716 


3018 


SC.path50 


2 


SO 


49 


190 


146 


190 


434 


434 


190 


190 


190 


190 


190 


190 


SC.pathSOO 


2 


SOO 


499 


1982 


1496 


1982 


7818 


7814 


1982 


1982 


1982 


1982 


1982 


1982 


SC.pathlOOO 


2 


1000 


999 


3980 


2996 


3980 


17706 


17726 


3980 


3980 


3980 


3980 


3980 


3980 


SC.path50 


7 


SO 


49 


114 


98 


114 


258 


180 


114 


114 


114 


114 


114 


114 


SC.pathSOO 


7 


SOO 


499 


1162 


998 


1162 


3754 


3260 


1162 


1162 


1162 


1162 


1162 


1162 


SC.pathlOOO 


7 


1000 


999 


2326 


1998 


2326 


7S62 


7114 


2326 


2326 


2326 


2326 


2326 


2326 


SCsimpleCyclcSO 


2 


SO 


SO 


202 


ISO 


202 


436 


452 


284 


284 


202 


202 


202 


202 


SCsimpleCyclcSOO 


2 


SOO 


SOO 


2000 


1500 


2000 


7782 


7804 


2968 


2968 


2000 


2000 


2000 


2000 


SCsimpleCyclelOOO 


2 


1000 


1000 


4000 


3000 


4000 


17742 


17746 


5964 


5964 


4000 


4000 


4000 


4000 


SCsimpleCyclcSO 


7 


SO 


SO 


120 


100 


120 


258 


178 


132 


132 


120 


120 


120 


120 


SCsimpleCyclcSOO 


7 


SOO 


SOO 


1170 


1000 


1170 


3770 


3256 


1328 


1328 


1170 


1170 


1170 


1170 


SCsimpleCyclelOOO 


7 


1000 


1000 


2334 


2000 


2334 


7S76 


7128 


26S6 


2656 


2334 


2334 


2334 


2334 



Table 3: Summary for the instances solved by a polynomial time algorithm. 
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known optimal solution 










Graph 


d 


n 


m 


OS 


DB 


NAM 


RAM 


RCAM 


G2 


BFSG 


TFSG 


CHLS 


PEHVNA 


SFHWI 


RGjrandomAS 


2 


SOO 


6126 


- 


48361 


98140 


97640 


97566 


92076 


96816 


98004 


86628 


89684 


93230 


RCrandomBS 


2 


SOO 


6175 


- 


48880 


99288 


98334 


98344 


92836 


97S56 


98610 


87540 


90S62 


94990 


RGjrandomAlS 


2 


SOO 


18654 


- 


201234 


299322 


298002 


297740 


290074 


297304 


298738 


281686 


286346 


292550 


RGj-andomBlS 


2 


SOO 


18887 


- 


204529 


302016 


301636 


301704 


293428 


300738 


302536 


285518 


290026 


295792 


RG.randomA2S 


2 


SOO 


31254 


- 


379064 


S00784 


499788 


499738 


490796 


499032 


500724 


480966 


486306 


49S964 


RG.randomB2S 


2 


SOO 


31114 


- 


376931 


497684 


497456 


497524 


487788 


496732 


498384 


478812 


483780 


493214 


RGj-andomA3S 


2 


SOO 


4360S 


- 


574180 


699352 


697740 


697594 


687438 


697148 


698420 


677870 


683514 


691738 


follow-up on the next page ... | 
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Graph 


d 


n 


m 


OS 


DB 


NAM 


RAM 


RCAM 


G2 


BFSG 


TFSG 


CHLS 


PEHVNA 


SFHWI 


RG_randomB35 


2 


500 


43595 


- 


574020 


699236 


697246 


697540 


686808 


696756 


698494 


677294 


683192 


690294 


RG_randomA45 


2 


500 


56653 


- 


782958 


908042 


907104 


906882 


896474 


906354 


907852 


886786 


892358 


899844 


RG.randomB45 


2 


500 


55627 


- 


766539 


891646 


890022 


890094 


879572 


889732 


891178 


870416 


876214 


884740 


RG_randomA55 


2 


500 


68320 


- 


978888 


1095540 


1093664 


1093718 


1083122 


1093128 


1094468 


1074734 


1079810 


1090610 


RGj-andomB55 


2 


500 


68701 


- 


985749 


1101882 


1100048 


1099958 


1089652 


1099576 


1101090 


1080722 


1085512 


1094074 


RG_randomA65 


2 


500 


81279 


- 


1212022 


1302766 


1301848 


1301356 


1291892 


1300882 


1302240 


1284086 


1288564 


1295348 


RGj-andomB65 


2 


500 


81172 


- 


1210096 


1301186 


1300226 


1299860 


1289980 


1299550 


1300660 


1282456 


1286966 


1293572 


RG_randomA75 


2 


500 


93347 


- 


1429246 


1496336 


1495498 


1495320 


1486398 


1495054 


1495980 


1479758 


1484062 


1490928 


RG_randomB75 


2 


500 


93399 


- 


1430182 


1497266 


1496448 


1495956 


1487202 


1496172 


1497070 


1480906 


1484748 


1492458 


RG.randomA85 


2 


500 


106047 


- 


1657846 


1699742 


1699524 


1699104 


1692306 


1698948 


1699578 


1687470 


1690196 


1693914 


RGj-andomB85 


2 


500 


106111 


- 


1658998 


1700626 


1700554 


1700062 


1692992 


1699822 


1700268 


1688546 


1691352 


1697034 


RGj-andomA95 


2 


500 


118499 


- 


1881982 


1899982 


1899538 


1899100 


1895412 


1899376 


1899696 


1893354 


1894696 


1897084 


RG.randomB95 


2 


500 


118606 


- 


1883908 


1900908 


1901264 


1900678 


1897246 


1900698 


1900804 


1895088 


1896222 


1898688 


RG_randomA5 


7 


500 


6126 


- 


21504 


40774 


46738 


40522 


38070 


39990 


40518 


36392 


37494 


39780 


RG_randomB5 


7 


500 


6175 


- 


21700 


41204 


47124 


40886 


38358 


40222 


40816 


36570 


37794 


40504 


RG_randomA15 


7 


500 


18654 


- 


84924 


124204 


142714 


123766 


119752 


123254 


124042 


117348 


118986 


124010 


RG.randomB15 


7 


500 


18887 


- 


86322 


125500 


144464 


125316 


121384 


124626 


125386 


118666 


120570 


125044 


RGj-andomA25 


7 


500 


31254 


- 


160524 


207686 


239158 


207686 


203046 


206856 


207712 


199806 


201802 


207458 


RG.randomB25 


7 


500 


31114 


- 


159684 


206620 


238134 


206588 


201778 


205850 


206552 


198944 


201050 


206506 


RG.randomA35 


7 


500 


43605 


- 


234630 


289994 


333908 


289584 


284438 


288968 


289716 


281458 


283704 


289994 


RGj-andomB35 


7 


500 


43595 


- 


234570 


290000 


333834 


289650 


284588 


288862 


289652 


281502 


283608 


290000 


RG_randomA45 


7 


500 


56653 


- 


312918 


376680 


433734 


376316 


371212 


375824 


376474 


368306 


370416 


376604 


RG.randomB45 


7 


500 


55627 


- 


306762 


369694 


426026 


369460 


364486 


368858 


369702 


361318 


363534 


369646 


RG_randomA55 


7 


500 


68320 


- 


382920 


454306 


523376 


454142 


448710 


453224 


453876 


446038 


448168 


454306 


RGj-andomB55 


7 


500 


68701 


- 


385206 


456812 


526276 


456536 


451350 


455894 


456454 


448314 


450618 


456812 


RG.randomA65 


7 


500 


81279 


- 


460795 


540234 


622664 


540098 


535146 


539482 


540146 


532718 


534564 


540234 


RG.randomB65 


7 


500 


81172 


- 


460159 


539620 


621798 


539254 


534380 


538760 


539478 


532046 


534046 


539620 


RG_randomA75 


7 


500 


93347 


- 


548776 


620442 


715076 


620622 


616240 


619842 


620570 


613984 


615588 


620442 


RGj-andomB75 


7 


500 


93399 


- 


549192 


621312 


715644 


620878 


616030 


620392 


620938 


614182 


616400 


621312 


RG.randomA85 


7 


500 


106047 


- 


650376 


705130 


812572 


705052 


701458 


704718 


704864 


699976 


701318 


705130 


RG.randomB85 


7 


500 


106111 


- 


650888 


705676 


812948 


705388 


701914 


705162 


705472 


700330 


701820 


705286 


RGj-andomA95 


7 


500 


118499 


- 


749992 


788218 


907730 


787964 


786042 


787888 


788060 


785158 


785984 


788214 


RG.randomB95 


7 


500 


118606 


- 


750848 


788782 


908742 


788638 


786674 


788496 


788652 


785996 


786804 


788642 


Pet03_randomAl 


2 


1000 


4974 


- 


29154 


89750 


88988 


88944 


80096 


87408 


86806 


71874 


75440 


86824 


Pet03.randomA2 


2 


1000 


24738 


- 


239917 


446298 


444500 


444384 


426856 


442944 


445890 


411242 


42540 


443786 


Pet03_randomA3 


2 


1000 


49820 


- 


577482 


897992 


895654 


895760 


873202 


894196 


897554 


852854 


864912 


894160 


Pet03_randomA4 


2 


1000 


8177 


- 


56759 


147424 


146646 


146528 


135366 


145032 


146066 


125374 


130530 


144756 


Pet03.randomG4 


2 


1000 


8173 


- 


56961 


147164 


146030 


146536 


98482 


93990 


96648 


74282 


106720 


135804 


Pet03.hcl0 


2 


1024 


5120 


- 


29696 


56320 


91468 


91672 


56320 


88684 


84266 


56320 


56320 


56320 


Pct03.mesh33x33 


2 


1089 


2112 


- 


8320 


18942 


41788 


38714 


19350 


26152 


23268 


18722 


18900 


18904 


Pct03.3clt 


2 


4720 


13722 


- 


63462 


169346 


328306 


313178 


189096 


174530 


220584 


120332 


132904 


168708 










fol 


ow-up on th 


c next page 















Table [4]- summary for the instances without a known optimal solution - follow up | 


Graph 


d 


n 


m 


OS 


DB 


NAM 


RAM 


RCAM 


G2 


BFSG 


TFSG 


CHLS 


PEHVNA 


SFHWI 


Pet03_airfoill 


2 


4253 


12289 


- 


56732 


148896 


293980 


273364 


165388 


153078 


188894 


106416 


119024 


148382 


PotOa.crack 


2 


10240 


30380 


- 


145618 


726664 


788288 


762606 


426592 


393516 


497116 




441042 


696122 


PetOS.whitakerS 


2 


9800 


28989 


- 


134741 


375730 


752132 


723426 


371946 


355240 


492390 


301320 


335630 


375298 


PetOS.big 


2 


15606 


45878 


- 


212875 


650814 


1190894 


1189886 


726976 


650746 


830690 


436908 


482936 


649748 


Pct03_wavc 


2 


156317 


1059331 


- 


6884189 


21067766 


36008138 


34921810 


23977688 


21016361 


23184961 






20711426 


Pet03.cly 


2 


828 


1749 


- 


8609 


24712 


31192 


30752 


21392 


26068 


22924 


16884 


19846 


23174 


Pet03.c2y 


2 


980 


2102 


- 


10246 


29726 


37394 


37392 


25282 


30232 


26722 


20478 


24110 


11592 


Pet03.c3y 


2 


1327 


2844 


- 


13578 


41996 


56492 


54426 


35066 


44664 


38692 


28810 


33736 


33736 


Pet03.c4y 


2 


1366 


2915 


- 


13529 


43490 


57848 


56106 


37034 


44186 


38306 


27930 


34124 


42814 


Pet03.c5y 


2 


1202 


2577 


- 


12120 


37636 


50712 


48414 


32894 


37942 


33524 


25572 


29328 


35858 


Pet03.gd95c 


2 


62 


144 


- 


643 


1016 


1384 


1354 


1080 


1024 


1002 


866 


916 


920 


Pet03.gd96a 


2 


1096 


1676 


- 


7021 


30310 


33124 


30774 


23050 


27908 


19060 


18004 


19926 


28526 


Pet03.gd96b 


2 


111 


193 


- 


971 


2410 


2246 


2200 


1762 


1820 


1768 


1486 


1760 


1576 


Pct03.gd96c 


2 


65 


125 


- 


495 


1276 


1394 


1236 


882 


936 


964 


824 


950 


844 


Pct03.gd96d 


2 


180 


228 


- 


1002 


2446 


3070 


2952 


2592 


2802 


2050 


1822 


2054 


2024 


Pct03_randomAl 


7 


1000 


4974 


- 


14006 


35988 


37952 


35680 


32204 


35032 


35120 


29574 


30608 


35002 


Pet03_randomA2 


7 


1000 


24738 


- 


96266 


178872 


189270 


178334 


171354 


177498 


178924 


166598 


169020 


178052 


Pct03.randomA3 


7 


1000 


49820 


- 


244920 


360072 


381548 


359372 


350714 


358698 


359846 


343664 


347592 


359368 


Pct03_randomA4 


7 


1000 


8177 


- 


26710 


59134 


62466 


58842 


54702 


58128 


58810 


51290 


52990 


58048 


Pct03.randomG4 


7 


1000 


8173 


- 


26697 


59124 


62296 


58840 


40294 


39452 


40402 


31838 


43626 


56006 


Pet03.hcl0 


7 


1024 


5120 


- 


14336 


26204 


39034 


36824 


26280 


35020 


34452 


25892 


25940 


26192 


Pet03.mesh33x33 


7 


1089 


2112 


- 


4224 


8294 


16044 


15298 


8326 


10478 


10076 


8224 


8262 


8286 


Pet03.bintreel0 


7 


1023 


1022 


- 


2044 


7384 


7752 


7288 


7384 


6488 


2856 


3030 


4418 


6436 


Pet03.3clt 


7 


4720 


13722 


- 


27694 


68636 


132340 


120844 


73754 


71042 


86290 


52420 


55328 


68576 


Pet03.airfoill 


7 


4253 


12289 


- 


24792 


60750 


118422 


107822 


64224 


63170 


78482 


46768 


49868 


60466 


Pet03.crack 


7 


10240 


30380 


- 


69741 


277928 


293238 


287434 


168802 


158460 


193856 


138178 


175252 


269060 


Pet03.whitakcr3 


7 


9800 


28989 


- 


58482 


154188 


279798 


273056 


153490 


147924 


191430 


126198 


137200 


153272 


Pet03.big 


7 


15606 


45878 


- 


92511 


261418 


442852 


442412 


286180 


264862 


317434 


186158 


199324 


260680 


Pct03.wavc 


( 


156317 


1059331 


- 


3299657 


8223278 


14474822 


13238898 


9203558 


8150052 


9050366 






8110512 


Pet03.cly 


7 


828 


1749 


- 


4075 


10224 


13294 


12326 


8778 


10426 


9454 


7508 


8290 


9652 


Pet03.c2y 


7 


980 


2102 


- 


4822 


12218 


16010 


15044 


10540 


12132 


11058 


8724 


9986 


11592 


Pet03.c3y 


7 


1327 


2844 


- 


6436 


16810 


21678 


20930 


14652 


17786 


15230 


12216 


13676 


16612 


Pct03.c4y 


7 


1366 


2915 


- 


6442 


17168 


22212 


21438 


14884 


17320 


15490 


12460 


13732 


16872 


Pct03.c5y 


7 


1202 


2577 


- 


5750 


15112 


19450 


18730 


13452 


15170 


13618 


10946 


11978 


14520 


Pct03.gd95c 


7 


62 


144 


- 


320 


474 


780 


608 


446 


460 


492 


410 


404 


450 


Pet03.gd96a 


7 


1096 


1676 


- 


3800 


12056 


12736 


12092 


9514 


11106 


7990 


7768 


8308 


11050 


Pet03.gd96b 


7 


111 


193 


- 


491 


1030 


1062 


928 


744 


716 


772 


682 


788 


696 


Pet03.gd96c 


7 


65 


125 


- 


250 


574 


680 


542 


426 


448 


448 


378 


438 


392 


Pct03.gd96d 


7 


180 


228 


- 


555 


1028 


1258 


1188 


1060 


1126 


876 


838 


886 


958 



Table 4: Summary for the instances without a known optimal solution. 



List of acronyms 

• OS = optimal solution (if known). 

• DB = degree bound. 

• NAM = normal arrangement. The vertices {vi, V2, . . . , v„} of the guest graph are mapped to the leaves of the d-regular tree in their canonical ordering, i.e. by 0(fi) = bi, for 
i = l,2,...,n. 

• RAM = random arrangement for k = 1000. k random mappings of the vertices of the guest graph into the leaves of the d-regular tree are constructed, their objective function 
values are computed, and the random mapping with the best objective function value is selected. 

• RCAM = random contiguous arrangement for k = 1000. k random contiguous mappings of the vertices of the guest graph into the leaves of the d-regular tree are constructed, 
their objective function values are computed, and the random mapping with the best objective function value is selected. 

• G2 = arrangement produced by the leaf-driven greedy heuristic, see Section |4. II 

• BFSG = arrangement produced by the breadth- first search based greedy heuristics which tries each vertex as the starting vertex, see Section 14.11 If the graph has more then 
one connected components, they are arranged in a random order. 

• TFSG = arrangement produced by the depth-first search based greedy heuristics which tries each vertex as the starting vertex, see Section l4.1I If the graph has more then one 
connected components, they are arranged in a random order. 

• CHLS = arrangement produced by the construction heuristic which uses the local search approach to solve the MCBSSP, see Section 14.21 

• PEHVNA = arrangement produced by the pair-exchange heuristic for vertices which starts with the normal arrangement, see Section 14.3.11 

• SFHWI = arrangement produced by the shift-flip heuristic which accepts non-improving shifts, see Section 14.3.21 The algorithm terminates if no improvement is reached after 
3 days of running time. 

• - = the solution could not be found in a reasonable amount of time. 



